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Abstract 

We extend the phase space slicing method of Giele, Glover and Kosower for per- 
forming next-to-leading order jet cross section calculations in two important ways: 
we show how to include fragmentation functions and how to include massive par- 
ticles. These extensions allow the application of this method not just to jet cross 
sections but also to cross sections in which a particular final state particle, includ- 
ing a D or i?-meson, is tagged. 



1. Introduction 



The identification ( "tagging" ) of one or more particular particles in the final state of 
reactions at colliders has been used succesfully in the past to enhance specific signals. 
A recent tremendous success was of course the top quark discovery || at the Tevatron, 
where b quark tagging was used to suppress the background. It is expected that future 
analyses will increasingly rely on tagging various kinds of particles, witness the large 
effort spent on developing Cherenkov ("RICH") and micro- vertex detectors. Heavy (D 
and B) mesons in particular will certainly continue to play an essential part in this. It is 
imperative that the theoretical side of such analyses keep pace with these developments. 
A key ingredient of this theoretical effort must be the construction of fully exclusive 
Monte Carlo type programs, which include exact higher order QCD corrections and 
appropriate fragmentation functions that parametrize the transition of a parton to the 
particular hadron being tagged. 

It is our goal to provide a systematic general-purpose method for constructing such 
Monte Carlo programs. In this paper we present the formalism. Numerical studies and 
applications will appear elsewhere f2j. 

Our method is an extended version of one already employed in the calculation of 
higher order jet cross sections, and is known as the "phase space slicing method" |3]||||. 
One alternative method for calculating general jet cross sections is the "subtraction 
method" |J. This method has also been applied to inclusive heavy quark cross sections 
0. Another, quite recent method is the "dipole method" ||, which can also be used 
for tagged reactions, and whose generalization to heavy quarks is in progress M. Yet 



another method, employing a small gauge boson mass, was used in [fLOH . The above 
methods differ in how they approximate the phase space and matrix elements in the 
neighborhood of divergent regions. 

The particular version of the phase space slicing method we extend is given in the 
papers of Giele, Glover and Kosower in Refs. and ||. Their version minimizes com- 
putational effort in the calculation of jet cross sections by using the concepts of color 
ordering and crossing. The latter property, implemented in higher orders with the help 
of so-called crossing functions, allows the calculation of the matrix elements to be per- 
formed with all partons in the final state. We will employ both concepts here, and 
introduce the final state counterparts to the crossing functions, "tagging functions" . In 
addition we compute the heavy quark contributions to the matrix elements, and to the 
crossing and tagging functions. The results presented in this paper have already partly 
been used in Ref. [^TJ, and were briefly presented in ]12| . 



The paper is organized as follows. In section 2 we review the basics of the phase 
space slicing method in the formalism of Giele, Glover and Kosower, and recall how it 
is used for computing general jet cross sections. In section 3 we extend their method to 
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include fragmentation functions. The extension to heavy quarks is given in section 4. 
We summarize and conclude in section 5. 

2. Phase Space Slicing Method for Jet Cross Sections 

In order to be self-contained and to establish notation used in later sections, we review 
in this section the phase space slicing (PSS) method, as given in Refs. for calcu- 

lating jet cross sections to next-to-leading order (NLO) in QCD. The method admits 
in principle an arbitrary number of jets in the final state, and is applicable to jet cross 
sections with either leptons or hadrons or both in the initial state. It allows for the study 
of any distribution of a given jet cross section, and for changes of jet-finding algorithms 
without analytical recalculation of matrix elements. 

Consider the situation where at lowest order one has n massless partons p in the 
final state produced by a vector boson V, i.e. the reaction: 

(ll'^)v^n V . (1) 

V itself can be the result of a lepton (Z) anti-lepton (J 1 ) collision. At this lowest order 
each parton is associated with its own jet, hence the matrix element for this reaction 
can be used to describe the process: 

(IV jets, (2) 

but also, after appropriate crossings and convolution with parton distribution functions, 
the reactions 

hh! -> V(-> 11') + (n - 2) jets , (3) 
Ih -> l'+ (n - l)jets, (4) 

where h and hi are hadrons. One may also consider reactions without the vector boson 

V in (|) and describe pure QCD reactions like: 

hh' -»• (n - 2) jets. (5) 

These reactions (0-^) represent the most prominent jet cross sections measured at col- 
liders. When one performs the phase space integral over the final state parton momenta 
numerically, it is clear that at lowest order, where each parton forms a separate jet, 
one may impose any jet-defining cuts or acceptance cuts on the calculation. Although 
jet-finding algorithms can also be implemented easily, such a calculation cannot repro- 
duce the details of the corresponding experimental cross section as it does not yield a 
prediction for the internal structure of a jet. 
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It is the main purpose of the PSS method as described in Refs. [§],[|] to maintain 
these features (i.e. the applicability of matrix elements to different processes, the easy 
numerical implementation of cuts and acceptance criteria, and the ability to change the 
jet algorithm without any analytical recalculations) at next-to-leading order. 

The reasons for incorporating next-to- leading corrections in the modelling of jet cross 
sections are well-known: 

• Reducing the normalization uncertainty via decreased sensitivity to changes in the 
renormalization scale and/or factorization scale. 

• Testing the convergence of the QCD perturbation series for the observable under 
investigation. 

• Modelling of inter-jet hadronic radiation. 

• Modelling sensitivity to the jet algorithm. 

The latter occurs because of the presence of the n + 1 parton matrix element in the 
next-to-leading corrections. It is now possible for two partons to form a jet. However it 
is up to the particular jet algorithm to decide whether, and how, to merge a given pair 
of partons into a jet. 

For definiteness we shall consider the process: 



where we assume that the vector boson has been produced in an e + e~ collision. For full 
details, we refer the reader to Ref. Q. Later in this section we describe how to make 
the transition to the other reactions ([HJ). That part of the discussion is a synopsis of 



At next-to-leading order, the contribution of the n + 1-parton matrix element to 
the n-jet cross section in e + e~ collisions involves integration over configurations where 
one parton is soft or two are collinear. In such situations divergences arise. They cancel 
against corresponding divergences in the virtual corrections to the n-parton contribution. 
The PSS method isolates these divergences by slicing the n + 1-parton phase space into 
"hard" and "soft plus collinear" regions. 

The hard region of phase space is essentially defined such that all of the n + 1 
partons are resolved: they are separated in phase space such that the invariant mass 
of any given two partons is larger than an arbitrary (small) theoretical cutoff s m i n , 
s ij = (Pi + Pj) 2 = 2Pi ' Pj > s min, where Pi and Pj represent the momenta of partons i 
and j. In this region of phase space the n + 1-parton cross section is finite and all phase 
space integrals may be done numerically. We may write this contribution as 



V — > qq + (n — 2)g 



(6) 



Ref. [§. 



da R (e + e 



(n + 1) partons) 



1 1 



M(l,...,n+1) dP R (Q;l,... 



n + 1) , (7) 



2 2 2Q 2 
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where the spin averages and the flux factor have been shown explicitly, and |M(..)| 2 
is the n + 1-parton LO matrix element squared and spin-summed. The resolved phase 
space in Eq. (^) is given by 

dP R (Q;l,...,n + l) = 1 f T[e(s ij -s aia )dP(Q;l,...,n + l), (8) 



n g m f n f q \Il f n f q \ 



'j 



where Q stands for the sum of the momenta of the incoming electron and positron and 
dP is the usual phase space measure in 4 dimensions: 

n+1 //3 n n+1 

dP(Q ; 1, . . . , n + 1) = [ n J^^j (2^(0 - £ Pi ) (9) 

An identical particle factor for gluons and (anti)quarks of the same flavor is included. 
The product over 9 functions is defined such that all non-divergent, finite contributions 
are properly included in this resolved contribution. See below Eq. ([[7]) for more detail. 

In order to understand the behavior of the matrix elements in the soft and collinear 
region it is very useful to decompose the amplitude for producing n + 1 partons in the 
final state using color-ordering |13|]. To illustrate this, let us consider here for definiteness 
the amplitude for n g + 1 = n — 1 gluons radiating off the colored quark line that couples 
to the vector boson V**, 

M(l, . . . , n + 1) C1C2 = V%{K\ 1, . . . , n g + 1; K) ClC2 , (10) 

where K, K are the momenta of the final state quark and antiquark respectively and 
Ci, C2 are their color indices. One may write 

S,(K; 1, . . . , n g + 1; K) ClC2 = ieg n * +1 £ (T ai . . . T^^S^K; 1, . . . , n g + 1; K) . 

P(l...n 9 +1) 

(11) 

Here S^K; 1, . . . , n g + 1; K) is a colorless subamplitude in which the n g + 1 gluons 
couple in a color ordered way to the quark line, and e, g are the QED and QCD coupling 
constants respectively. The full amplitude is the sum over all permutations of these 
gluons. 

A similar decomposition may be performed for the amplitude with four or more 
quarks in the final state. We shall not discuss such amplitudes in any detail. We merely 
note that such amplitudes are color-suppressed in relation to the gluonic ones, and can 
be straightforwardly included in the formalism (see Ref. || for details). 

Using Eq. (pi]), the square of the full amplitude can be expanded in the number of 
colors N. One finds, after contraction with the vector boson V, and factoring out an 
overall factor (N 2 -1)/N 



V»S,. 



e 2 (V) (^ v -)(^% + 1 ] + ^[% + l] + ---), (12) 



2 J v N J\ 1 y 1 N 2 
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where 

A[n g + l}= £ \V»S ll (K;l...n g + l;K)\ 2 . (13) 

P(l...n s +1) 

There is no known general expression for the term B[n g + 1], but for any given n g it 
may be computed straightforwardly |TJ]]. Note also that Eq. ( 12 ) is a finite series. For 
example, with n g — 1 only the A [2] and £?[2] terms are present, and 5 [2] reads 

B[2] = -\W( S,(K; 1, 2; 7T) + 2, 1; K) )| 2 . (14) 

Here the 3-gluon coupling cancels, and the gluons behave in fact like photons [14 . 



We now examine the complementary soft plus collinear region where one or two 
of the Sij are smaller than the theoretical cut-off, s mm . In this region only n of the 
7i+l partons are resolved. The divergences that occur in this region are regularized by 
working in d = 4 — 2e dimensions. The advantage of the color-ordered decomposition 
of the matrix element is that an n + 1-parton colorless ordered subamplitude factorizes 
into an n-parton amplitude and a 1-parton factor, for an appropiately small value of 
s m i n . The n + 1-parton phase space also factorizes into the n-parton phase space and 
a 1-parton factor. The integration of the 1-parton amplitude factor over the 1-parton 
phase space can easily be performed analytically. 

After this integration the squared soft and collinear matrix elements, resp. (V^S^H 
and \V»S, 



(they are given explicitly in Eqs. 



25]) and (pq)), can be combined with the 
leading order (LO) and NLO virtual (V) n-parton contributions to form an effective n 
resolved parton squared matrix element: 

|M(1, . . . , n)| e 2 ff = |M(1, . . . , n)| 2 G+v + \V%? 8 + \V»S,\ 2 C . (15) 

After coupling constant renormalization |M(1, . . . , n)| 2 ff is finite, by virtue of the renor- 
malizability of QCD |15| and the KLN theorem Jl6| . 



The resolved n parton contribution is now given by: 

1 1 



da R (e + e 



Note that because 



sure here in 4 dimensions. 



n partons) 



2 2 2Q : 



M(l, 



n) eS dP H (Q;l, 



,n) 



(16) 



M(l, . . . , n) is finite, we may take the n-parton phase space mea- 



Combining Eqs. flTB] ) and (I7|), we may now write the next-to-leading order n-jet cross 
section in e + e~ collisions as 



da(e + e — > n jets) = 6 



da R (e + e 



n partons) 



da R (e + e 



[n 



1) partons) 



(17) 

Here G represents all the experimental effects, including the jet algorithm. The integral 
in the second term represents the projection of the phase space of the n + 1 resolved 
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partons onto the n jets phase space. We can now be more explicit about the definition 
of the product over the ^-functions in Eq. (|8|), which is implicit in the second term 
in Eq. (O). For each color order contributing to |M(1, . . . , n + 1)| 2 , the product runs 
over no more than two invariants involving the same momentum, otherwise some finite 
contributions of order s m i n / s m i n might not be properly included. 

Eq. ([H]) now involves only finite quantities and may be used to construct a fully 
differential NLO Monte Carlo program for the production of n jets in e + e~ collisions. 
Note that both terms on the right hand side of Eq. (|TTD depend on s min , but that 
no observable may depend on it, as long taken small enough such that the 

approximations used in the soft and collinear region are valid. For a given process this 
must be checked numerically with the Monte Carlo program. 

We now briefly describe the procedure to calculate the soft and collinear contribu- 
tions. 



2.1. Soft Behavior of Matrix Element, Phase Space, and Cross 
Section 

Colorless ordered subamplitudes share the property of QED amplitudes that, in the limit 
where one of the external gauge bosons becomes soft, the amplitude factorizes into a 
lower order amplitude and an eikonal factor: 

S^K; 1, ...,n g ,s;K) -> e x (n g ; s; K) S^K; 1, n g ; K) , 
S^K; 1, ...,m, s,m + 1, ...,n g ;K) -> e x (m; s; m + 1) S^K; 1, n g ; K) , 

S„(K;s,l,...,n g ;K) -> e\(K; s; 1) S^K; 1, n g ; ~K) , (18) 

where the eikonal factor is given by 

e A (a;s;6)= . (19) 

V Pa ■ Ps PfPs J 

Here e^(s) is the polarization vector of the soft gluon with momentum s and helicity A, 
and the momenta p a and pb are those of the two hard partons adjacent to the soft gluon 
in the specific color order considered. Summing over the helicities of the soft gluon one 
obtains for the square of the amplitude 

I^mI 2 - ^(^n^T 1 )! E s F (K- 1 l,...,n g ;K)\V^(K;l,...,n g ;K)\ 

P(l...n 9 ) 



2 



2 / V N 

i_ 

' N 2 ' 



+ 0(±)], (20) 



where 

<g 2 N 



s F (K; l,...,n g ;K) = \e x (K; s; 1)| 2 + . . . + £ \e x (n 9 ; s; K)\ 2 ] . (21) 

A A 
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In the soft limit the <i-dimensional (d = 4 — 2e) n + 1-parton phase space measure, 



n+l rid— I 

dP d (Q;l,...,n+l) = []\ 



U2^1 (2 ^ (Q -? P - ; 



also factorizes, up to terms of order s min /Q 2 , according to 



dP d (Q ; 1, . . . , a, s, b, . . . , n + 1) -> dP d (Q ; 1, . . . , n) dP^ oft (a, s, b) 



(22) 



(23) 



where 



dPsoft(a,s,b) 



(4*0 



i67r 2 r(i 



Safe 



SasSbs 



Sab 



°as min Sbs) • (24) 



The ^-functions define the soft region formally. After combining (]20D , (21) and (24), 
and integrating over s as and s& s the effective squared matrix element with n g resolved 
gluons and two resolved quarks is 



AT 



F(l...n s ) 



iV 2 ' 



n,;ir)|^(K;l,...,n 9 ;K)| 2 

(25) 



with 



5 F (^;l,...,n 9 ;/f) = (— — A—) +--- + (7 i 



(26) 



Here /i is the dimension regularization scale, and a s = a s (fi) = g 2 fi~ 2e /47r. We may 
interpret the dimensional regularization scale /i here already as the renormalization scale 
fin, as coupling constant renormalization is implicitly assumed. Note that the identical 
particle factor in Eq. (§), l/(n s + l)!, is changed to 1/%! in Eq. (p3|), thereby taking into 
account the (n g + 1) possibilities of choosing a soft gluon. The contribution to | V i S' At || 
at higher order in 1/N can be calculated straightforwardly. 



2.2. Collinear Behavior of Matrix Element, Phase Space and 
Cross Section 

The matrix elements also develop a singularity when two of the partons, say a and b, 
become collinear, i.e. 

Pc=Pa + Pb] Pa = Zp c ; p b = (1 - Z)p c 

Sab =PI= Zpa-Pb < S min - (27) 
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where s m i n is again assumed to be small. In this case the full squared matrix element 
factorizes according to [] 



1^(1, 



,a,b, 



,n + 1; 

It is useful to extract color factors via 



^(l,...,c,...,n) 



xc—*ab 



9 2 N 



r 



Then we have 



c-^ab 



P. 



c— >a&\ 



Sab 



(28) 



(29) 



(30) 



The splitting functions P^ a b rea d, in the conventional scheme (in which both external 
and internal particles are treated in d = 4 — 2e dimensions) 

/l + z 4 + (l-z) 4 \ 



P. 



g^qq 



P, 



q->qg 



(z) 



V z(l-z) 
2 f z i + (l- z Y 
N\ 1-e 



2(1- 



* 2 -e(l-z) 2 



(31) 



Notice that only when two adjacent partons in the subamplitude in a specific color order 
become collinear does one find a divergence. Here, the four quark (+(n g — 1) gluons) 
amplitude, T^, must be considered: when a neighboring quark and antiquark form a 
flavor singlet and become collinear, then this amplitude also contributes to the n parton 
amplitude. In analogy to Eq. (|20"D, we have therefore in the collinear limit: 



V»T U 



2( g 2 N Yg( N*-U 



P(l...n g ) 



N 

,n g -K)\V»S»{K-l, 



where 

c F (K; 1, 



K^qg 



+ f 



f K ^ 9q + n g n f f 



,n g ;K)\ 2 + 0( 



N 2 ' 



g~*qq 



(32) 



(33) 



The last term in this expression is due to 7) t . The factor n g in front of this last term 
arises because the identical particle factor in Eq. (H) for the four quarks contribution 
changes from l/{n g — 1)! to 1/%!. The factor rif results from summing over all (light) 
quark flavors. In the collinear region, defined in Eq. (p7[), the phase space measure also 
factorizes, up to terms of order s min / Q 2 as: 



dP d (Q ; 1, . . . , n - 1, a, b) -> dP d (Q ; 1, . . . , n - 1, c) dP c e on (a, b; z) 



(34) 



1 In Ref. || the notation c^ b ^ c was used to indicate the splitting of parton c, instead of c c f* ab . 
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where 



(47T) 



-ds a bdz 



s ab z(l - z) 



v^mm Sab) 



(35) 



167r 2 r(l - e) 

To obtain the collinear behavior of the cross section one needs to integrate the following 
expression: 



9 2 N 



f^ ab dPUa,b;z) 



a s N\ (Airfi 



2\e 



2tt mi - e 



d>S a b S ab 



1 f 1 -** 



4, 2l 



dz 



z(l-z) 



p 



c^ab \ 







>r(i-«) 


■'mm 



_ Ic-*ab\ z li z 2) 

e 



(36) 



where the boundary values Z\ and z 2 prevent the integral from picking up contributions 
from the soft region: 



S a -lo > z l s a-lc 



s bb+l > (1 — ^2) ^cb+l — Smin ■ 



(37) 



These boundaries clearly depend on the color order of the particular subamplitude 
under consideration. The collinear behavior of the effective squared matrix element with 
n resolved partons, of which one is an unresolved collinear pair, is now given by 



V»S„ 



v 2 1 v n 

x [ £ C F (K;l,...,n 9 ;K)\V^S,(K;l,...,n g ;K)\ 2 + 0( 

P(l...n g ) 



N 2 ' 



(38) 



where 

C F (K;l,...,n g ;7C) 



C( S N\ 1 /47T/i 2 Nel 

~2Wrn -e) vZ~T^ 7 



(0, 



(39) 



j- / ^min "Smin \ , 



ski S12 



Note that when the boundary corresponds to a soft fermion, there is no singularity, so 
we can put the corresponding Z{ to zero. For completeness we list here the results for 
the I functions, defined in Eq. (^), for massless partons J|||. 



(zi,z 2 ) 



zi +z; 



11 .7T 2 67n 2 

h ( e + 0{e 

6 V 3 18' V 



Iq-* qg ( z l, Z 2 ) = (1 - J^)[( 



Z2 6 ~ 1 



3 ,7T 2 7x 

4 + ( ¥-l) e 



0(6 2 



(^1,^2) 



1 rl 5e 



iV 



3 + 9 



+ 0(e 2 ). 



(40) 

(41) 
(42) 
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For I q ^ qg (zi, z 2 ) the 1/N 2 term in fact contributes to the 0(1/N 2 ) terms in Eq. (j3S|). We 
use this definition of I q ^ qg because it is more convenient for what follows. The remaining 
J's are trivially derived from these, via the relations 

Iq-+gq\ z ll Z V ~ Iq-^qg\ z 2-> z l) > Ig-*qq[ z l, z 2) = Ig-^,qq\ z 2i z l) (43) 

and 

Iq->qg( z l, z 2) = Iq-*qg\ z l, z 2) , Iq-*gq\ z li z 2) = Iq-*gq( z l, z 2) ■ (44) 

Again, we refer to Ref. [Q for all details. 



2.3. Crossing 

In this subsection we review how to describe reactions involving initial state hadrons 
with the PSS method. The original derivation and all details can be found in Ref. ||. 

The calculation of the NLO matrix elements in the previous section was done with 
all partons in the final state. To describe processes with initial state partons at NLO, 
one would like to generalize the crossing property of the LO matrix elements to NLO. 
In particular, one would like to use the already calculated effective matrix elements 
of Eq. (pUD , and do the crossing numerically, without having to redo the analytical 
calculation for the specific process under consideration. This is possible, via universal 
"crossing functions" ||. The feasibility of thus extending crossing to NLO rests again 
upon the simultaneous factorization of phase space and matrix element for both initial 
and final state collinear radiation. 

In general the NLO differential cross section for a process with initial hadrons % \ 
and "K 2 may be written as 

d^Ki5c a = J2j dx iJ dx 2 3f 1 (x l )3* i \x 2 )d<jl LO (x l ,x 2 ) , (45) 

a,b 

where a, b denote parton flavors and x\,x 2 are parton momentum fractions. The sym- 
bols in this equation do not have quite the same meaning as they do in the standard 
factorization formula [17] . Here 9^ (x) is an "effective" parton distribution function and 



da^ h LO is not the next-to- leading partonic cross section (which we denote by da^ b LO ). 
Rather, it consists of the finite effective all-partons-in-the-final-state matrix elements, 
in which partons a and b have simply been crossed to the initial state, i.e. in which 
their momenta p a and Pb have been replaced by —p a and — p&. The contributions that 
distinguish da^ b LO from da^ b LO are in fact included in S^fix) and are what distinguish 
the latter from the usual parton distribution functions f^(x). Explicitly, we have 

do N ah LO = da L ° + a s d(6a» LO ) + 0(a% (46) 
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with a s = a s (fj,R) from the crossed matrix elements, and 

= fFfaup) + a.C?(x,p F ) + 0(a 2 s ), (47) 

where mass factorization has taken place at scale /j, f . f^(x,^ F ) is a NLO parton 
distribution function and C"^(x,ijlf) is a finite, process-independent crossing function. 

There are two contributions that distinguish da^ b LO from the full next-to-leading 
partonic cross section da^ b LO . First, when crossing a parton into the initial state we 
must correct for the possibility that this parton is in fact an unresolved collinear pair. 
Such "fusing" processes are not to be taken into account at leading twist and must hence 
be subtracted. Second, the contribution from the initial state collinear region must be 
included. Schematically, one has therefore, before mass factorization, 

-El/' -/c M (-)W*) - f?V) C dzP a ^ b {z)] S ^ , (48) 
^ l Jx z z Jo J e 

where the /^ are the bare distribution functions and the collinear singularities are still 
present. After mass factorization in a particular scheme, e.g. MS, one may write 

C?' m (x, n F ) = f [A? (x, M ln(^E) + B?** (x, „ F )] . (49) 
Ztt l fi z F J 

The functions A and B can be expressed as convolution integrals over the parton distri- 
bution functions and are listed in Ref. p. The B functions are scheme dependent. The 
crossing functions, being process independent (universal), need to be calculated only 
once for any set of parton distribution functions (see [|I1||T9| , [20| for recent sets). 

The full NLO differential cross section can now be written as: 
dax 1 n 2 = J^Jdxtj dx 2 f^ 1 (x 1 ,iJ F )f^ 2 (x2,fXF)da a N b LO (x 1 ,X2) + 

a,b 

+ a a (fi F ) (C^ 1 (x l ,i2 F )f b ^ 2 (x 2 , fj, F ) + f^ 1 (x u /2 F )C b K2 (x 2 , ii F ))d<j a L b °(x 1 ,x 2 ) 



a,b 



where the scheme dependence of the individual functions is suppressed. 

We have reviewed quite extensively the PSS method for the calculation of general 
jet cross sections involving massless partons, and now proceed to the main task of this 
paper, which is to generalize the method to include fragmentation functions (in the next 
section), and heavy quarks (in section 4). 



3. Fragmentation and Tagging Functions 

The PSS method can be extended to describe not just jet cross sections, but also 
cross sections in which a particular particle is identified or "tagged" in the final state. 
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Here that particle is assumed to be a hadron K . Descriptions for such processes require 
the introduction of a fragmentation function D^f(z), representing the probability of a 
parton h to fragment into K , with z the momentum fraction of h carried by K . In this 
section we show how to include fragmentation functions in the PSS method. 

Let us consider the NLO cross section for the production of a hadron K and (n — 1) 
jets in e + e~ collisions. 



daK=Y. d °h L0 ^(z)dz, (51) 
h 

where Dj^ is an effective next-to-leading order fragmentation function, in a similar sense 
that 3^ in Eq. (|45|) was an effective next-to-leading order parton distribution function. 

is in fact given by Eq. ([17]), with one of the resolved partons fixed to have flavor 
h. The sum is over all parton flavors h in the intermediate state. The goal is again to 
use the results already obtained for the production of n jets at NLO without need to 
modify the effective matrix element for the resolved n-parton contribution (Eq. ([IB]) ). 
We may expand 

dal LO = dal° + a s d(5a% LO ) + 0(a 2 s ) , (52) 
<Df (z) = D« (z) + a s T? (z) + 0(a 2 s ) , (53) 

where the are the final state equivalents of the crossing functions, described in the 
previous section. We call them "tagging functions" . 

Substitution of Eqs. ( [52] ) and ( |53"D into (|5T| ) gives up to order a s 

da^ = £ [{da L h ° + a s d(8a» LO ) )ti* (z) + a s da L h °T^ (z)} dz . (54) 

h 

What follows is in fact very similar to the derivation of the crossing functions in Ref . || . 
The tagging functions consist of those corrections that render Eq. (|54]) correct. 

In the evaluation of the contribution in the collinear region in section 2, the inte- 
gration over the momentum fraction z was performed. This collinear integration should 
however have been done in convolution with the fragmentation functions. Therefore, 
the integrated contribution must first be subtracted, and then the properly convoluted 
contributions should be added. For a specific color ordering, the integrated contribution 
is given by: 

(*) = £/ ^WpUp, «; y)Dh (*) , (55) 



p,u 



so that 



i3* (*) = -(St^(^) £ ;E Wi/i.w.)^W. (56) 

where the sum is over all allowed parton flavors. The Ih^ pu are given in Eqs. (|40[ - [44[) . 
Here y\ and 7/2 are the equivalent of the color-order dependent boundaries Z\ and z 2 in 
Eq. (|37|). The complete tagging functions T^ 1 will be independent of the color ordering. 
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Next, we need to add the collinear terms having the proper convolution with frag- 
mentation functions. For each specific color ordering, the resulting coefficient is 

« s Tf c „ m {z)dz = £ / dPUp, «; y)#^(y)i>? (x)6{z - xy)dx , (57) 

p,u J 

where y = E p /Eh and x = E<k / E p , so that z = Eh /Eh is an energy fraction. Using the 
expressions of section 2, we find 

(58) 

Other definitions of z are equivalent, e.g. using the plus components rather than the 
energies. We have 

N 1 /47T/Ael ^ fl-Va dyv v1 -el 



(59) 

This may be written as 

W = -(l^^Ej; d f<^(y,v^(,M , (60) 



with 



i^(y, y 2 ) = ~y- £ (l - y)- £ P,%Jy)0(l - y - y 2 ) . (61) 

The functions Kh^ p (y, y 2 ) are equal to the functions J p ->h(y, y 2 ) in Ref. || in the crossing 
case, but multiplied by y~ e and with h and p inverted. Using the plus prescripton 

F(z)+ = lim (6(1 - z - p)F(z) - 6(1 - z - (3) dyF(y)) , (62) 

0->O v Jo ' 

we may expand the Kh^ p (y, y 2 ) in order to extend the y- integral to 1. For example, the 
right hand side in Eq. ( |60|) contains terms that may be expanded as 

Now combining the integrated and convoluted components leads to 
Th ( z ) = T h conv (z) — T h int (z) 

= E t —Y h ^ p (y)D M p (z/y) . (64) 

P Jz y 
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The Yh-> p (y) are independent of the color ordering-dependent boundaries y\ and 2/2, 
so that Tf^ is indeed an overall factor multiplying the LO partonic cross section, as 
advertised in Eq. (|54|). 

The functions (z) in (E4T) still contain collinear singularities. These singularities 
are cancelled via mass factorization applied to the fragmentation functions, which in the 
present case amounts to writing 



nCK .scheme / 



„ Jz y 



(65) 



where is the final state factorization scale. The Rh^ p (y, Hd) are counterterm func- 
tions, the finite terms of which determine the factorization scheme. In this section, 
where we are only dealing with massless partons, we choose the MS counterterm func- 
tions Rh^p, as in Ref. M for the crossing functions: 



N\ (Ann 2 



1 Uri 



2tt/ V fijj J r(l-e) e2 



1 . 



1 + y 2 
1-3/ 



N_\ 

27 ) l7!T, 



1 1 



.Z2pt = (y) 



2^ 



1 1 



V? D J r(l-e)e4 



g^qq 



(y) 



(66) 

(67) 
(68) 



2tt, 



i4> 



1 l[ (lliV-2n / ) 



+ 2 



l-z) 



+ 



r(l-e)r 

(1-*) 



6N 



8(1 -y) 



+ z(l-z) 



(69) 



where P^ib can be inferred from the massless splitting functions in Eq. ([Tl]). The 
counterterm functions in Eqs. (p^-p9]) cancel the 1/e poles and the factor (47r) e /T(l — e), 
present in Y^ p , and replace the dimensional regularization scale [i by the final state 
factorization scale fin- Then we have 



1 h \ z -> 



^D) = Y,f-{ Y ^y) + Rf%(y^D)]Df(z/y^ D ) (70) 
t) y v / 



Notice that the effective fragmentation function D is not dependent on the factorization 
scheme or scale 



(z) = D% (z, fi D ) + a s T? (z, p D ) + 0(a 2 ) . 
The results of the calculation in Eq. (|70| ) are summarized by 



2vr 



U£{z^ D )H^) + V^ u \z^ D ) 



(71) 



(72) 
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where 



U? (z, Ato) = £ (*, Hd) , Vf MS (z, p D ) = £ V^ S (z, M . 



(73) 



The functions U%Z p are scheme independent and are given by: 



Ug^g (Z,H D ) 



* y 



(UN-2n f ) 



+2/(1-2/) 



QN 



D?(z/y,n D ) 



5(l-y) + 2 



(1-2/) 



(1-2/)- 



N 2 Jz y 

1 Z" 1 ^/ 

* 2/ 



3 AM \ i 1+2/ 2 

2 5(1 ^ ) + (T^); 



D Q K (z/y, f M D ), (75) 



U* g (z, nn) = \t -P e q Z" 9q {y)Df (z/y, p D ) , 



(74) 



(76) 



(77) 



where the functions PyZ^ pu are obtained from Eq. ([31]) . The scheme-dependent functions 



y^H,MS can ^ e s j m ji ar jy derived: 



9^3 1 



2,/i D ) 



1% 

* 2/ 



2 In 2/1 



(1-2/)- 



+ ^ + 2/(1-2// 



/V 67 5n/\ r/1 . n An(l-y) 

+ + — 1 \S(l-y) + 2y(— — 

V 3 18 9NJ v y; y V i _ 

+ 2 (2-=^+y(l-y)) 111(1-2/) 



(z/y, fi D ) 



(7S 



V; 



s 1/ 
= o(l 



Of, MS 

q->q v 2 V ' iV 2 



i r 1 ^ 

* 2/ 



hi 2/ 



1+2/ 2 



,7T 



7. 



(l-2/) + + ( T-2^-^ 



+ (l-y) + (l + y 2 ) 



2x,ln(l-2/) 



1-2/ 



D?(z/y,pL D ) 



V g ^(z, M = \f z < ^ j - \p;Z° q M In (2/(1 - 2/)) " iW2/) 



q ^ S (z^n) = U 1 - 



D?(z/y,fM D ) 



4J Z y 



P q Z° gq (y) in (2/(1 - 2/)) - ^(y) 



Df(z/y,ii D ). 



The functions -P/_^ jfc are defined by 



PI 



PiZ" jk + ePL jk + 0(e' 



(79) 

(80) 
(81) 

(82) 
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with the Pf_+j k listed in Eq. fl3ip. Notice the In?/ terms generated by the extra factor y~ e 
in the collinear phase space measure. We stress that the tagging functions are universal 
and need only to be calculated once for a given set of fragmentation functions (see |21| 
for recent sets). 

We can combine the crossing functions and tagging functions in order to describe a 
NLO fully differential cross section for collisions of hadrons "K \ and "K 2 in which hadron 
"K is tagged in the final state: 

damx^n+x = E / dx 1 dx 2 dz^ 1 (x 1 )^ b H2 (x 2 )da^°(x l ,x 2 ,z)^(z) . (83) 

a,b,h 

Substituting 

fodbh°( l ?i,x 2 , z) = daa° h (x 1 ,X2i z) + a s d5a^° (x u x 2 ,z) + 0(a 2 ) , (84) 

3V (*) = f? (x, » F ) + a .C% (x, Hf) + 0(a 2 s ) , (85) 
(z) = (z, /in) + a s T™ (z, y. D ) + 0(a 2 s ) , (86) 



leads to 

da-K^^ji+x = 2J 



^MfFfa) (da™(x 1 ,x 2 ,z) + aJ8a» L h °(x 1 ,x 2 ,zj) D%{z) 



abh 

+a s C^ l {x l )f h Ji2 {x 2 )da L ab ° h {x ll x 2) z)D^ (z) 
+a s f^ 1 (x 1 )C^ 2 (x 2 )d^ b ° h (x 1 ,x 2 , z)D% (z) 

+a.f™(x 1 )fj*\z a )do*°{z 1 ,z 2 , z)T? (z) + 0(a 2 ) 



dx\dx 2 dz , (87) 



where all the matrix elements in the partonic cross sections (da^ b % da^°) are crossed 
versions of the ones with all resolved partons in the final state. 

This completes the treatment of the massless fragmentation functions. We next 
extend the PSS formalism to include heavy quarks. 



4. Massive Particles 

A very important class of tags consists of heavy (D or B) mesons, indicating that a 
heavy quark was produced in the hard scattering. A description of such reactions with 
massless quarks is often not sufficiently accurate. In this section, we show how to extend 
the previous results to cases involving heavy quarks with an explicit mass. Throughout 
this section the mass is understood to be renormalized in the on-mass-shell scheme. 
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4.1. Soft Contribution 



We consider the reaction: 



(«'-)v-QiQ 2 + K + l)0, 



(88) 



where V stands again for any vector boson (e.g. 7, Z, W, etc). Now Qi and Q2 are 
heavy quarks with masses m\ and m 2 , respectively. Our results can be straightforwardly 
extended to the case with four massive quarks V — > Q1Q2Q3Q4 + (% — 1)<?, or two heavy 
and two light ones, etc. 

When gluon s become soft, the colorless subamplitude again factorizes into an 
eikonal factor multiplying the ordered subamplitude with one less gluon. The square of 
the eikonal factor is, summed over the helicities of the soft gluon, 



fab(s) = ^|e A (a;s;6)| 



2 4s ab 4m 2 4m l 



Sas^bs 



2 ' 



where we keep the definition of s as as in section 2: 



(89) 



(90) 



and m a and m& are the masses of the neighboring particles, which can be mi, m 2 , or 
zero. To analyse the soft behavior of phase space, we only need to understand the three 
particle case, without loss of generality. The three particle phase space measure reads 
in d = 4 — 2e dimensions 



. , , W W/'' V, A .,, 

2E a 2Eb 2E. 



dP d (Q;a,b, s) = (2tt) 



S d (Q -Pa-Pb- Ps) 



(91) 



The soft region is defined as in the massless case: s as ,Sb s < s min , and s min is assumed 
to be much smaller than Q 2 . This phase space measure again factorizes, up to terms of 
order s min /Q 2 : 

dP d (Q ; a, b, s) - dP d (Q ; a, b) dP^a, b, s) . (92) 
The soft gluon measure 



dP soft( a ' b > S ) = ( 27r ) 1 ~ d of 8 ^( s min - Sas)0(s mi n - S bs ) 



2E„ 



(93) 



may be written as 



dPsoft(a,b,s) = 



(4vr) e A (e " D 



22 22 

Sas Sfe S a fe — m b S ab ~~ m a S bs 



167T 2 T(l - e) 
ds as dSi, s 9(s m - m S as )$(s m ; n Sfe s ) , 



X 



with 



\ 2 /i 2 2 
A = s ab -4m a m b . 



(94) 
(95) 
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Note that if one takes m a and m b to zero, one recovers the massless result in Eq. (p4|). 
The factor between square brackets is in fact proportional to (pjr) 2 , where pj: is the 
component of p s in the direction perpendicular to p a in the center of mass frame of a, b, 
and s, and must hence be positive. This positivity condition leads to the kinematical 
constraint Q 

r_s bs < s as < T+s bs , (96) 

or equivalent ly 

a_s as < s bs < a+s bs , (97) 

where 



Sab 

T± = 7T-?.± 



Sab \ 9 ^o,b , 

) 9 > a ± = o o ± 



2m 2 \ 2m 2 m b ' 2m^ ^ 2m 2 a 

Note that (t_|-) _1 = a T . For the rest of this discussion we assume m a > m b . 

One can show that r_ < (>) 1 if s a t, > (<) m 2 a + m 2 . Furthermore r + > 1. The 

complete integration region for s as and is conveniently divided into two regions. The 
first region is 

s bs E [0, cr_ 

■Smin] j S a s e [T-S bs , T + S bs ] , (99) 

and the second region is 

1) s ab > m a + m b : s bs E [o"-S m i n , s m ; n ] , s as E [T-S bs , s m - m ] , 

2) s ab <m 2 a + m 2 b : s bs E [o"_s min , cr+s min ] , s as E [r_s bs , s min ] . (100) 

The threshold condition is s ab > 2m a m b . Note that m b = implies cr_ = 0, in which 
case the first region of integration does not contribute, and in the second region the 
lower limit of s bs reduces to zero. 

Combining the soft part of the phase space and matrix element in this general massive 
case, we find 

S(m aj m b ) = ]{m U s)^Ms). (101) 

We distinguish between two cases in the rest of this section, because an extra divergence 
occurs when m b = 0. We find the following results after performing the integration in 

Eq. grog): 

m b = case 

Denoting m a = m we find 

S(m ,0) = a -f 1 [*OLX (^)'j im ,0) (102) 

27T r i-£ S min / VSmin/ 



2 This result can also be derived from the condition that | cos(6>)| < 1, where 9 is the angle between 
a and s in the center of mass frame of a, b, and s. 
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where, corresponding to the two ranges in ( |100| ) 
1) s ab > m 2 : J(m, 0) 



2e 2 \m 2 



2) s ai , < m 



J(m,0) 



Sab 



2e \m 2 
1 1 
2^2 + 2e 



1 



C(2) 



Note that the strength of the 1/e 2 poles is the same in both cases. 



2 

Sab 



(103) 

(104) 



m b 7^ case 

Here we find 

S(m a ,m b ) 



'47T/i 2 



2\l-2e 



2tt r(l-e) 



A2" e 



J(m a ,m b ). 



(105) 



The contribution to J coming from the first region of integration (|99"D has the following 
form: 



J(m a ,m b ) 







2e 






-T-) 







2(r_ - r+) + (r_ + r+) ln(^) - 2(r+ - r_) ln(r + - r_) 



^±^ln(^)[l + 21n(r + -r_)] 



Li 2 (i_z±)_Li 2 (i-— ; 



(106) 



The contribution from the second region of integration (ITDO) is 

1) s afe > m 2 + m\ : J = 1 - (r_ + r+) + r_r + 

+ ln(r+)(r_ - r+ - (r_ + r+) ln(r_)) + |(r_ + r+) m 2 (r + ) 

2) ^<m 2 +m 2 : J = (r_ — r + ) ln(— ) + il—hZ±l m 2 (— 



(107) 
(108) 



We cannot take the limit m b — > in this expression, which reflects the extra singu- 
larity that leads to a double pole in Eq. ( |104|) . Note finally that if m a = m b = m there 
is only the case s ab > 2m 2 (and a± = r±), because the other case violates the threshold 
condition. 



It is obvious how to modify the soft factor Sp in Eq. (|26|) for reaction 



n g > 



x 



n g = 



S F (Qi, l,...,n g ; Q 2 







V 2tt > 


r(i- e ) ( 



mm 

1 (S(n B -l)l 



(109) 



x 



°mm c °mm c "mm 

(m 2 ,) 1 2e ( m\ \e i 

/A Jmin 



+ J(0,m 2 )(^) 



(110) 
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Note the replacement of N by (N 2 — 1)/N in the case n g = 0. 

There are no 1/e singularities associated with the collinear region when one of the 
two collinear particles is a heavy quark. The mass of the quark screens the singularity. 
Therefore the above extension of the soft iT-factor to the massive case is all that is 
needed to be able to apply the PSS method to processes involving heavy quarks. As far 
as UV renormalization is concerned in this situation, any scheme is of course allowed, 
but the Collins- Wilczek-Zee [2^] scheme, in which the decoupling of the heavy quarks 
for small external momenta is manifest, is preferred. In this case only the light quarks 
contribute to the running of the QCD coupling. 

The PSS method, extended as above to include soft radiation from heavy quarks, has 
in the recent past already been used in constructing NLO programs for reactions where 
heavy quarks are produced in leptonic collisions [2^] , and also in hadronic collisions using 
crossing functions p4j|. We shall next extend it further to include collinear radiation. 



4.2. Collinear Behavior and Massless Limit 

In a scattering process involving a heavy quark in which the typical hard scale Q (e.g. 
the pt of the heavy quark) is of the order of the heavy quark mass m, the results 
obtained thus far are adequate to describe the reaction. Singularities due to collinear 
configurations involving the heavy quark are screened by its mass, logarithms ln(Q /m) 
are small, and the cross section is well-behaved. However if Q/m is large, the collinear 
region contributes large logarithms In (Q/m) that dominate the cross section at higher 
orders in perturbation theory. Moreover, they may spoil the convergence of the QCD 
perturbation series and one may wish to resum them. The ln(m) terms are the equivalent 
of the 1/e poles that appear in the collinear region with massless partons in dimensional 
regularization and can be treated conceptually on equal footing. In what follows we 
describe how heavy quarks can be incorporated into the PSS formalism if Q ^ in. 

The outline of this subsection is as follows. We first present results for the behavior of 
the cross section in the collinear region when the heavy quark is in the final state in the 
limit that m/Q and m? / s m \ n are small. Then, we derive the heavy quark contribution 
to the tagging functions. We present results for two factorization schemes, a minimal 
substraction ("MS") and massive MS("MMS") scheme, the latter defined such that 
the heavy quark fragmentation functions are the same as in the MS massless case. We 
show that the counterterm functions of the massive MS scheme are in fact equal to the 



perturbative fragmentation functions of Ref. [25]. The fragmentation functions absorb in 



this way, through factorization, the ln(m) terms and resum the ln(Q /m) terms through 
their evolution. We then present the equivalent results for the initial state. 

A comprehensive treatment of heavy quark cross sections should be performed in the 
context of a variable flavor number scheme f2~6|| , which aims to describe the full region 
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from Q ~ m to Q ^> m. This subject has recently been an active topic of investigation 
27| , |28| , p9| , |30| , |3Tl , |32| . The results presented in this paper can be used to implement this 
scheme and we comment on this possibility at the end of this section. 



Final State 



We first derive the behavior of the cross section in the collinear limit when the typical 
hard scale Q is much larger than m. As in the massless case, we take all partons and 
massive quarks as outgoing. Because the mass serves as a regulator for the collinear 
divergences, we may work in four rather than in d dimensions. We follow the structure 
of section 2.2. Let us denote the particles that become collinear by A and B, and their 
parent particle by C, 



Pc 
Sab 



Pa+Pb ; Pa = zpc ; Pb -- 
™c = (Pa + Pb) 2 ~ m 2 c < s r 



;i - z)pc 



P 2 c 



(111) 



5ab is the inverse propagator of C. The capital letters for the momenta remind us that 
two of the A, B and C particles are heavy quarks. We have defined z in Eq. (|1 1 1|) as a 
4-momentum fraction. Other definitions (such as an energy, or 3-momentum fraction) 
are equivalent in the collinear region when m < Q. As usual we assume that s min is 
small compared to Q 2 . We furthermore neglect terms that give contributions of the 
order of m 2 / 's min to the cross section. 



As in the massless case, in the collinear region the matrix elements factorize: 



1^(1, 



,A,B, 



n + 1) 



c 



C^AB 



,C, 



n 



with 



and 



C^AB 



f 



c 



C^AB 



C^AB 



C^AB 



z,m,5 A B) 



6 



(112) 

(113) 
(114) 



AB 



The massive splitting functions Pq_^ ab contain a term involving the mass m of the heavy 
quark (Q) and Sab- 



P^ Qg (z,m,5 Qg ) =2(1 



1 . 
iV 2 "' 



1 + 2? 



2m 2 



PQ^ gQ (z,m,5 Qg ] 
where 5 Qg = 2p Q ■ p g (= s Qg ), and 



/ ! + (! 



z 



5 Qg 
2m 2 



8, 



Qg 



2 

N 



z 2 + {\- zf + 



2m 2 



QQ 



(115) 
(116) 

(117) 
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where 5qq = (pg + pq) 2 - The terms containing the mass can be calculated by the 
methods used in |33| . 

The phase space also factorizes in the collinear region: 

dP\Q ; 1, . . . , n - 1, A, B) -> dP A (Q ; 1, . . . , n - 1, C) dP coll (A, B; z) , (118) 

where 

dP coll (A, B; z) = ——dzdb AB . (119) 

l07T^ 

The integrated collinear factor is then given by: 

J(^f)f C - AB dP coll (AB;z) = (^-)Ic^ab(z u z 2 ), (120) 

where 

I C ->ab(zi, z 2 ) = - dz -ABpg^fa m , Sab) (121) 

4 J 2i J<5 min OAS 

Zi and z 2 are defined as in the massless case, and are needed to avoid the soft region 
when A or B is a gluon. The kinematic limits 5 m i n are different for each of the splittings 
and are given in Table [I]. Note that although m 2 is small compared to s m j n , we must 



Table 1: The kinematic limits S m \ n for tagging functions. 



C -> AB 


^min 


Q^Qg 


m a (l— z) 

z 


Q^gQ 


rri^z 
(l-z) 


g^QQ 


m 
z(l-z) 



retain it in <5 min , because it will give rise to In (m 2 / s m ; n ) contributions to the cross section 
that cannot be neglected. Other mass terms in the limits of integration that give rise to 
order (m 2 /s min ) contributions have been neglected. We define at this point the K c ^a 
functions via: 

Ic^ab{z\,z 2 ) = dzK C ->A(z,z 2 ). (122) 



When z 2 cannot be put to one (when B is a gluon) its effect is absorbed in the integrand. 
With the plus prescription defined in Eq. (|62]), we obtain the following results: 



1 1 

K Q ^ Q (z,z 2 ) = 



z 2 )ln[ S -^z\ -2z) (123) 



1 — z) + ^ v ' \ m 2 
(1 + - 2 )( ! 4^) + ~ '(» - «)(21n*(ln ] - 1) - ln 2 « 
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l + (l-z)\ f s min l-z^ 2(1 -z) 



ln( 



m A z 



v ' ; 2N 



(z 2 + (1 - zf) In ^fz(l -z))+ 2z(l - z) 



(124) 
(125) 



Consequently, 



(0,0) 



f dz 

2N Jo 



(z 2 + (1 - z) 2 ) In ( - z)) + 2z(l - z) 

In ( ^ m ' n ^ 
2N V3 9 



(126) 



1 1 r 1 -^ i + -- / 



1 ~t~ 2! N , /-S m j n ^ \ 2 - (127) 



m 2 \ — z> 1 — z 



1/ 1 X 



-21nz 2 --)ln(^) +ln 2 z 2 + 21nz 2 - — + - 



and 



(Z!,0) = Jq^q^O,^). (12? 



The massive Ic^ab functions can be used, as in the massless case, to calculate the 
collinear factor Cp{. . .) (see Eq. (|3DD) for each of the colorless subamplitudes. For suf- 



ficiently inclusive, infrared-safe observables, the ln(m 2 ) terms of the Ic^ab functions 
will cancel with corresponding terms of the virtual contribution, the same way the 1/e 
terms cancel in the massless case. In these cases it is better to integrate the collinear re- 
gion analytically as done here and perform the cancellations analytically. The numerical 
integration and cancellation, although always possible, is in practice difficult at large Q . 

We now turn to the problem of including the heavy quark into our formalism for 
fragmentation functions. For this we must add the heavy quark contribution to the 
tagging functions, where we assume that the transition of the heavy quark to a heavy 
hadron, at least for values of Q much larger than m, may be described by convoluting 
the partonic matrix element with a (non-perturbative) fragmentation function. In the 
normalization conventions of this section the heavy quark contribution to the tagging 
functions can be constructed via 

(z) = (z) - Tl mt (z) = £ f ^-Y H ^ P (y, m)D^(z/y) , (129) 



where 



T^iz) = (^Z/^W^)^-), (130) 
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and 



1 HAnt 



(131) 



J{can now be a light or heavy hadron and H, P, and U can be heavy quarks. 
There are four possible contributions due to the heavy quark: 



•2tt' 



K, 



Q-*Q 



(y, 2/2) - 5(1 - y)I Q ^ Qg (0, y 2 ) 



v 2ttV iV 2; 



1 + 2/ 2 
1-2/ 



In 



+ 



;i + y 2 )ln 2/ 



Yg^ Q (y,m) 
Y Q ^ g {y,m) 



(i-y)+ 



ln(l - y) 



,7T 



(1-2/) 
2 3 



Wl- 2/) 



(2/) 



Z7T 



and 



iV 



„(y, m) = -(—)I g _> QQ (0, 0)5(1 - y) . 



(132) 

(133) 

(134) 
(135) 

(136) 



Yg^g is the heavy quark contribution that must be added to the massless Y g ^ g . As in the 
massless case, all dependence on the boundary y 2 has cancelled. Eq. Q53|) can now also 
be applied to the massive fermion case. The heavy quark contribution to the tagging 
functions is independent of the hard process. 

When ln(Q /m) becomes too large we need to these logarithms to all orders. This 
can be done by the evolution of the fragmentation functions, through factorization. The 
factorization procedure, as in Eq. flTCP, leads to scheme dependent Y H ^ P functions: 



Y H ^ P (y, m) -> y~(2/, fi D ) = Y H ^ P (y, m) + ir~(y, ™, fi D ) . 
A minimal substraction, "MS" , choice of the counterterm functions Rh^p is: 



pMS 



. ,N.l. 1 , 



( 1±2A 



1 ( m2 \ 

H'D 



N 1 m 2 



L Q^g\ 

R™ Q (y,m,ii D ) 



'-)-p^°-( v )H— 



N 9 m 2 

RT^ s (y,m,^ D ) = _(_)_in(^. )S(1 - y) 

Zir bJM [i D 



g^g 



(137) 

(138) 

(139) 
(140) 
(141) 
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where the P^ab can be inferred from the massless splitting functions in Eq. fl3~ip . R™_, g 
is the massive contribution that must be added to the massless R g ^ g - Notice that the 
heavy flavor is active in the evolution of the fragmentation function. The Yjf^ p are now 
finite for m — > at fixed 

We can decompose the heavy quark contribution to the tagging functions as in 
Eq. (j72|) for the massless case, with: 



U, 



'K 



1 

2< 



J_. 

N 2 ' 





\C +y2 )] 


L y 





Dq{z/v^d) 



4 Jz y 



u, 



m,1K . 
9^9 



y 

= 1 f l dy_ 

iJ z y 

Z,fX D ) = - 
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~6N 9 



(142) 

(143) 
(144) 
(145) 
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V'q (z,h d ) 



and 



1 

2iV 



2/ 

J_, 

jV 2 ' 



y 2 + (l-y) 2 ) \n(y(l-y)) + 2y(l-y) 



D%(z/y,fi D ), (147) 



V, 



y 



l + {l-yf 



ln(- 



<A 2(1 -y) 



9^9 



D*(z/y,ti D ). 

(148) 
(149) 



The U are scheme independent, and have the same functional dependence on their 
respective fragmentation functions as in the massless case. The contribution of the heavy 
quark to U g -> 9 and V g ^ g increases the number of flavors n/ in the massless expressions 
by one. After factorization the heavy quark contribution to the tagging functions is 
independent of m and is therefore finite as m goes to zero. 

However, the tagging functions, although now m- independent, are not the same as 
those corresponding to the massless case of section 3. This implies that in the MS scheme 
the m — > limit of the full partonic cross section, which can safely be taken, is not the 
same as if the calculation had been done with massless partons, and in the massless 
MS scheme, from the outset. When Q ^> m, the most appropriate factorization scheme 
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is therefore one where the subtracted tagging functions are the same as in the MS mass- 
less scheme, up to small terms of order m 2 /Q 2 . In other words, the counterterms in 
this massive MS ( "MMS " ) scheme must be such that Ypf^p (z, /ip) is equal to the corre- 
sponding massless function Y^ p (z, fin)- From this requirement we can straightforwardly 
infer from our earlier results what these MMS counterterm functions must be, and we 
find 

AM , 1 



1 +V 2 \ , ,m 2 1 + y 

Mtt) + 



W=£^) + <l + »VM<l-») 



l-Z/ y+ 3d (1-3/) + 

(150) 




R^ S 9 (y,m,^ D ) = a^d/)|ln [ ^f- ] +1 ] , (151) 



R^ (y, m, fi D ) = R?* Q (y, m, fi D ) (152) 



R^ MS (y, m, p D ) = R^ s (y, m, p D ) (153) 

The choice of scheme has in turn implications for the fragmentation functions, which 
absorb these counterterm functions. They must be determined in other reactions in the 
same scheme. Notice from the above explicit expressions for the counterterms that the 
heavy quark contribution to the kernels for the evolution equations of the fragmentation 
functions are still those of the massless splitting functions. 

We remark that when the contribution from the collinear region is calculated analyt- 
ically in the MMS scheme as described in this section, the massive counterterm functions 
are actually not needed. All that is needed to actually perform a calculation, besides 
the MMS fragmentation functions, are the tagging functions, and in the MMS scheme 
they are known by definition from the massless MS case. However, when the contribution 
from the collinear region and the convolution with the MMS fragmentation functions are 
computed numerically, then it is necessary to explicitly add the counterterm functions 
in Eqs. (|150[ - |153| ) to make the calculation consistent. In that case the equivalent tagging 
functions that should be used in Eq. fl54|) are given by: 



£ f ^R^ s P (y, m, Vd)D¥ (z/y, M . (154) 
p Jz y 

The original approach towards resumming the large logarithms in heavy quark pro- 
duction at high Q 2 is due to Mele and Nason in Ref . [^5] , and is called the heavy 
quark perturbative fragmentation functions (PFF) approach. Here the large logarithms 
ln(Q jm) are summed at the partonic level. The resummation again occurs via evolution 
of the following initial distributions at scale /i (of order m) to the large scale Q 

2tt 



D Q (z,to,m) = 5(1 -z) + ?^d$(z,to,m), (155) 
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and 



D g (z,Ho,m) = ^ d { g ; (^,/i ,m) 



(156) 



In our formalism the d$ are actually precisely given by the counterterm functions (mul- 
tiplied by — 27r) given in Eqs. (|150| - |153|) , at [ir> — /-to- The expressions in Eqs. ( |150| ) and 



|152|) may be rewritten to obtain the same result as in Ref. p5 



N 2 -l 
2N 



(157) 



z 2 + (l- zf 



ln(4) 
rrr 



(158) 



The PFF approach has been used to obtain resummed transverse momentum spectra 
for a variety of reactions p^j35 . 



Initial State 

In this section we present the results needed to include processes with heavy quarks in the 
initial state within the PSS formalism, when Q > m. As in the massless case we assume 
that the calculation of the squared matrix element was done with all partons, including 
the heavy quarks, in the final state. The heavy quark contribution from the collinear 
regions is computed analytically, as just described, under the assumption that Q is large. 
We now calculate the heavy quark contribution to the crossing functions. We follow the 
methods and notation of Ref. ||, as outlined in section 2.3. Recall that we must add the 
contributions from the initial state collinear region and subtract the contributions from 
the final state collinear region that are part of the analytical calculation but that do not 
contribute when the unresolved collinear pairs in the colorless amplitudes are crossed. 
As in the final state case, we neglect terms of order m 2 /s min . 

The process under consideration may be generically represented as 

P + a^U + M, (159) 

where P and a are the initial particles of the hard subprocess and M can decay into 
any number of particles. Particle U is radiated collinearly off particle P (and is not 
observed) : 

P^H + U, (160) 

whereas H undergoes a hard scattering with a to produce the final state labeled M. 
Here, up to two of P, H, and U may be heavy quarks, hence the upper case notation. 
The momentum fraction z is again defined by: pn = ^VpiVu — (1~ z )Vpi see Eq. ( |1 1 1| ) 
and the comments below it. The inverse propagator of the off-shell particle H is given 
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by: 5 



PU 



(Pi 



PUj 



m 2 H . The collinear region is defined by \5pu\ < s m[n . In this 



limit the phase space and the flux factor factorize in the following way: 

—dP{P + a^U + M)^ dP coll {P, U) x -^—dP{H + a -> M) . (161) 

Note that sua = zsp a , where the z factor is included in dP co u(P, U). In our notation f], 
sp a = 2p a .pp = (p a + pp) 2 — nip = A 1 / 2 ((p a + pp) 2 , mp, 0), where we have assumed that 
m a = 0. 



The collinear phase space measure is given by: 

1 



dP coll (P,U) 



167P 



z dzd\5 



PU\ 



(162) 



The full squared matrix element factorizes in the collinear region as described by 
Eq. ( p8[) in the massless case. The dp~^ HU are now given by: 

V 9. >z 



(163) 



2 >z \5 PU \ 

This result can be obtained by appropriate crossing of the final state cp-factor in 
Eq. ( |113p . The massive splitting functions are given in Eqs. ( |115 - 117|) . The result 
for c^ QQ m Eq. (|163|) is also given in Ref. |29| . 

The initial state component of the crossing functions is given by 



a s { ~'H,init 
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1 dz 



x Z 



fp(x/z)J P ^ H (z, z 2 ) 



(164) 



Here "K is an initial state hadron, fp~ (x) is the unrenormalized distribution function of 
P in Ji . Z2 is introduced as before to avoid the soft singularity when U is a gluon. The 
kinematic limits 5 m \ n are given in Table ^. 

Using Eq. (|62|) , we obtain 

5(1 -z)(21nz 2 (ln(- Sr 
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(165) 



X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 
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Table 2: The kinematic limits S m i n for the crossing functions. 



p — > hu 


^min 


Q^Qg 


m 2 (l — z) 


Q^gQ 


a-z) 


g^QQ 


m 1 

a-z) 



Note that the Jp^u can also be obtained from the final state Kc->a by appropriate 
crossings. 

The final state contributions to the crossing functions are given by 



N 



H ,final{x) = ( — )fH(x)Y, I B^Pu(z 1 ,Z 2 ) 

Z7T „ 



(166) 



where the Ih^pu are given in Eq. ( |126[ - |128| ). These contributions have to be subtracted: 



( X ) — ^H,init ( X ) Ch, final ( X ) • 

Defining the functions Xp-^H(z,m) by 

C£ (x) = £ t -fp (x/z)X P ^ H (z, m) , 

r, J X Z 



(167) 



(168) 



the four heavy quark contributions are given by, 



X Q ^ Q (z,m) 
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(-)-(l--^ 



X g ^ Q (z,m) 
X Q ^ g (z,m) 
X™(z,m) 



(— ) In(^) 

* _(l+,>)(!=&^>) +( 

(1 — z) + v 1 — Z ' + 

N 

Z7T 

AT 

-^^qq(o,o)5(i-^). 



7T 



;)*(!-*) 



(169) 



where X™_^ g is the massive quark contribution to the massless X g ^ g . 

After factorization, using the minimal subtraction scheme defined by the massive 
counterterm functions given in Eq. (|138[ - [L41j) (and replacing \i D by fip), the crossing 
function can be written as in the massless case: 



C 



IK, MS, 
H 1 



x, m, fip) 
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(x, fi F ) ln( : ^) + B H ' M (x, p, F ) 
fi F 



(170) 



29 



with 



we have 



1 r 1 dz 

A^ g (x^ F ) = -J -fg(x/z^ F )P*Z° 9q (z) , 



1 r 1 dz 



A^ g {x^ F ) = -—f^{x^ F ) 



(171) 

(172) 

(173) 
(174) 
(175) 



and for the scheme dependent functions 

Bq^q(x^f) = ^(1 - Jfi) f ^fq{x/z,n F ) 
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B ^Q S ( x ^f) = ^ ^ ^-f^(x/z,fi F ) 



[z 2 + (1 - z) 2 ) ln(l - 2) + 2z(l - z) 



, (177) 
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X, *f(a;,// F ) = ^(1 - — ) j( ^-f£(x/z,n F ) 



l + (l-z) 2 ln( l~^ _ 2(1 -z) 



5 



78) 
(179) 



Just as for the tagging functions the massive quark contributions to A g ^ g and B g ^ g are 
equivalent to a change of by one. 



It is easy to verify that the choice of the massive MMS counterterms in Eq. ( |1 50j - |l 53| ) 
yields crossing functions that are equal to those of the massless case (given in Eq. (3.39) 
and (3.40) Q) in Ref. |. 

In both factorization schemes the large logarithms ln(Q /m) due to initial state 
collinear radiation are resummed via the evolution of the parton distribution functions. 



Matching at low Q 

The PSS method, with the extensions presented in this paper, provides a natural frame- 
work for the implementation of a variable flavor number scheme (VFNS) [ P^ , P^ , p8] , pP|j3l]| , 



4 We note a typo in Eq. (3.40), the sign after the delta function in the expression for B^j^ 18 should 
be changed. 
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3~T| , |3"2| in the calculation of differential cross sections. The name refers to the fact that in 



such a scheme the number of active flavors depends on the value of Q . Here we merely 
sketch the outlines of such an implementation, we hope to provide a complete example 
in the future 0. For a review of recent VFNS developments see e.g. Refs. pl| , |36[| . 

When the typical scale Q is close to the mass of the heavy quark, we may use 
the PSS method extended by the soft gluon radiation results of section 4.1. There is no 
need to compute analytically the contribution of (approximately) collinear configurations 
involving the heavy quark, a numerical integration should suffice. The heavy quark is 
not treated as a parton, neither in the initial nor in the final state, i.e. it does not 
contribute to the running of a s , nor to the evolution of the parton distribution functions 
or fragmentation functions. To calculate the production cross section of a D or B meson, 
we have to convolute the parton-level results with a non-evolving fragmentation function 
(e.g. a Peterson- Zerwas function, or simply D^(z) = 5(1 — z)). This calculational 



scheme, when adopted irrespective of Q , is called a fixed flavor number scheme (FFNS). 

When Q is much larger than m, the fully extended PSS method should be used, 
including the heavy quark crossing and tagging functions, factorized e.g in the MMS 
scheme, so that large logarithms ln(Q /m) are resummed, see earlier in this section. Now 
the heavy quark is treated as a parton, and it contributes to the evolution of a s , the 
parton distribution functions and fragmentation functions. 

The issue is now how we should merge these two descriptions into one calculational 
scheme. For the discussion here we assume that the factorization theorem holds even at 



a low scale close to the mass of the heavy quark, see Ref. [31]. We take /zr = //p = /i 
for simplicity. At the end we shall put \i — Q . 

Let us concentrate on the initial state. We have in mind a cross section sufficiently 
inclusive to not require any fragmentation and tagging functions^. The corresponding 
matrix element is computed with all partons in the final state, including the heavy 



quarks. Renormalization is performed in the CWZ scheme [22|, and all soft and collinear 



singularities have cancelled due to the KLN ]16[ theorem. 

We start by choosing a matching scale, /io — m, such that, for fi < /i , we use rii 
active parton flavors in our calculation and treat the heavy quark as non-partonic, while 
for fi > fi we use n\ + 1 active flavors by treating also the heavy quark as partonic in the 
sense that it participates in the evolution of a s and the parton distribution functions. 

For fi < /io we cross, following Eq. fl50|), all possible (rii) light partons to the initial 
state, and fold the result of each crossing with the corresponding n^-flavor parton dis- 
tributions functions. We also add the contributions from the corresponding (MS ) light 
parton crossing functions. The evolution of the parton distribution functions involves 
only the light degrees of freedom. 



3 E.g. one might think of the deep-inelastic heavy-quark structure function. 
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For fi > hq we replace all n^-flavor parton distribution functions and a s by their 
ni + 1 versions. The relations between these n\ and ni + 1 quantities are given by so- 
called matching conditions, computed for a s to two-loop in Ref. and for the parton 



distribution functions to one and two-loop in Ref. p9| , p7| . We extend the sums over 
parton flavors and their crossings in Eq. (JS^) to include the heavy quark. However, where 
we should add a heavy quark crossing function, we only need to add the corresponding 
counterterm function (e.g. those in Eqs. (|150| - |153| )). This is because, in the matrix 



element, all contributions, including near-collinear ones, are computed numerically []. 
Of course, for /i ^> ^ one may use the heavy quark crossing functions derived earlier 
in this section to compute the contributions of the collinear region and counterterm 
functions. 

To show that this calculational scheme constitutes a VFNS, we must verify that both 
the low and high Q descriptions are properly reproduced. Clearly, when /i ^> /i , we 
have an «j + 1 flavor description, where the large logarithm ln(Q /m), converted by the 
counterterm functions into ln(Q //z), may be eliminated from the matrix element and 
resummed into the parton distribution function by the choice fi = Q . 

For /i just above /io the requirement is that the cross section computed in the VFNS 
is the same as computed in the FFNS[]. The argument essentially follows Ref. [26|| . For 
fi > /io, the difference with the FFNS cross section is mainly due to the evolution of 
the gluon distribution function caused by gluon splitting into a heavy quark pair, and 
the counterterm functions. For [i ^ m (putting /i = m, and suppressing all irrelevant 
arguments), the former contribution can be written as 

cto-^CVFNS) - AxSi(FFNS) ~ (a s In (£) ~P S _^ da^ (180) 



2tt/ \mj 2 



with fi not much larger than m. The corresponding counterterm functions (see e.g. 
Eq. ( |140| )), precisely cancels the 0(a s ) term in Eq. (|180|) , as required. 



5. Conclusions 

In this paper we extended the phase space slicing method of Giele, Glover and 
Kosower for computing NLO corrections to jet cross sections to incorporate fragmen- 
tation functions and heavy quarks. This makes the method applicable to any reaction 
in which a particular final state hadron (light or heavy) is tagged. 

6 For \i = Q not too much larger than /io, neglecting terms of order m 2 /s m i n when computing heavy 
quark crossing and tagging functions, is not very accurate, due to the constraint s m { n «9 2 . Here it is 
therefore better to compute contributions from the collinear region numerically. 

7 Note that we consider one loop matching here, in which a s and the parton distribution functions 
are continuous across /io, in the MS scheme. 
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The extension to fragmentation functions for massless partons amounted to a gener- 
alization of the crossing function approach || to the final state, and led us to introduce 
a new class of functions we named tagging functions. 

The extension to heavy quarks consisted of two parts. Effects of soft radiation off 
heavy quarks were straightforwardly included. Contributions from collinear radiation 
involving heavy quarks were included in heavy quark crossing and tagging functions. We 
showed that heavy quark tagging functions led naturally to the Mele-Nason perturbative 
fragmentation functions [[25]. We described briefly how our results might be included in 



a variable flavor number scheme. 

In this paper we developed the formalism of the extended PSS method. Its practical 
use, and its application in the context of a variable flavor number scheme, will be assessed 
in future work 0. 
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